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1. Introduction 

For many years a wide variety of physical systems have been described, often 
approximately, in terms of coupled two-state systems P3 12 E]. In more recent years 
application has been found in quantum information and quantum computing 0j, where 
such two-state systems have been used to describe a quantum mechanical version of 
the classical computer bit. This quantum bit, or 'qubit', is described as a linear 
superposition of two states (say 'off' and 'on'), so that before a measurement the qubit 
is in some sense simultaneously both off and on, unlike a classical bit which is always 
either off or on. 

While two-state quantum systems are widely used, their utility suffers because 
there are only a limited number of known analytic solutions. For most two-state 
systems numerical calculations are required. Although Standard numerical methods 
are both fast and reliable for these simple systems, finding how the corresponding 
physical systems work is largely a numerical fishing trip in cloudy wàters. Analytic 
solutions, where they exist, are more transparent. One such solution that has received 
widespread applications in quantum òptics |U IHj is obtained using the rotating wave 
approximation (RWA). In this approach, periòdic transfer of the population within 
the system is achieved by applying an external field (e.g. làser) that is tuned to a 
narrow band of freqüències to match a particular transition between the system's levels. 
The result is well known Rabi oscillations. In this paper we wish to call attention to 
another analytic solution for two-state quantum equations, namely the limit of sudden 
pulses. Such a fast pulse is called a 'kick'. Unlike a periòdic field with well defined 
carrier frequency used in the RWA technique, a kick is localized in time and consists 
of a broad range of frequency components. Due to their 'non-periodic' structure, kicks 
are well suited for applications that require occasional modifications of the system's 
state (one kick - one transition). There are advantages in using kicks in systems where 
two states of interest lie close to one another in energy (nearly degenerate systems). 
And, of course, kicks would be a natural choice in cases where a pulsed source has 
to be used for one reason or another. Fast pulses are an essential ingredient in the 
kicked rotor or Standard map, a paradigm of the transition to chaos in one-dimensional 
time-dependent dynamics [6 . The kicked rotor was first realized in the laboratory by 
exposing ultracold sodium àtoms to a periòdic sequence of sharp pulses of near-resonant 
light 7 . Signatures of classical and quantum chaotic behavior, including momentum 
diffusion, dynamical localization, and quantum resonances have all been observed in 
such atom òptics experiments. Intriguing connections have also been demonstrated 
between momentum localization in the quantum kicked rotor and Anderson localization 
in disordered lattices [H] . The use of short pulses for the purpose of control of quantum 
systems was suggested previously for a variety of systems, including excitation of 
electrònic states in molècules jü] , product formation in chemical reactions [Tül HI] , and 
quantum computing [Ï2*] . Due to complexity of those systems, control pulses have to 
be carefully shaped to achieve effective control, and determination of such shapes often 
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requires one to use either numerical techniques, or genètic algorithms. The response 
of quantum systems to fast pulses has also been studied extensively in the context of 
pulsed nuclear magnètic resonance (NMR) [T3j. In this case, the pulse width is typically 
short compared with the time scales of the relaxation processes, including spin-lattice 
relaxation time T\ and transverse relaxation time T 2 * associated with line broadening. A 
single pulse may be used to study free induction decay; more sophisticated multi-pulse 
sequences are used in spin echo experiments and in multi-dimensional fourier transform 
NMR ^3]. Techniques similar to ones discussed here are therefore developed to study 
the detailed evolution of a single spin or múltiple spins, corresponding to multi-qubit 
systems, including spin precession and relaxation effects in the time intervals between 
pulses. Quantum gates necessary for computation have been constructed using such 
pulse sequences (TH], and realized in experiments [TH]. In some systems, discussed in 
this paper, one may take advantage of single or múltiple external pulses of simple shape 
(e.g. Gaussian) that contain many frequency components to reliably control quantum 
systems. There are other limits in which analytic solutions may be obtained, including 
perturbative and constant external interactions, and degenerate systems. However, 
these tend to be of limited use, as we discuss below. 

In this paper we develop simple analytic solutions for singly and multiply kicked 
two-state quantum systems |17j . Since the kicked limit is an ideal limit of pulses very 
sharp in time, we do numerical calculations for pulses of finite width to illustrate the 
region of validity for fast pulses. We do this for 2s — 2p dipole transitions in hydrogen 
and illustrate the limits on the band width of the signals required to sensibly access such 
a region. Part of our motivation for this study is to understand reaction dynamics [TH] 
and coherent control ^Hj in the time domain. In particular we are interested in the study 
of observable effects of time ordering and also in understanding the related problem of 
time correlation (2Ü| in few body dynamics |21j . corresponding to a system of a few 
dynamically coupled qubits. We use both analytic and numerical solutions to study 
effects of time ordering in multiply kicked systems. In a simple analytic example we 
illustrate the difference between time ordering and time reversal invariance. Atòmic 
units are used throughout the paper. 

2. Dynamics of a two-state system 

2.1. Basic equations 

A two-state quantum system may be described by a wave function, = 

° 2 í 1 J ' w ^ ere ^ x ) an< ^ ( ) re P resen t the two basis states, e.g. on and off. Here a\ 
and Ü2 are complex probability amplitudes restricted by the normalization condition, 
|ai| 2 + | a-2 1 2 = 1. 

The basis states I j and í j are eigenstates of an unperturbed Hamiltonian, 
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H , given here by, 

H = ~<T S , (1) 

where AE = E 2 — E% is the energy difference of the eigenstates of Ho, and a z is a 
Pauli spin matrix. The average energy of 'on' and 'off' states of the unperturbed system 
may always be taken as zero since a shift in overall energy of the system corresponds 
to an unphysical overall phase in the wavefunction. Probability amplitudes evolve as 
a,j(t) = üj(0)e~ lEjt , and the occupation probabilities of the basis states, Pj = \dj\ 2 , 
remain constant in time. 

The states of a qubit can be coupled via an external interaction V(t), so that the 
occupation probabilities change in time. For simplicity, we assume that the interaction 
has the following form: 

V(t) = V(t)a x , (2) 

i.e. all of the time dependence in the interaction operator V(t) is contained in a single 
real function of t and the interaction does not contain a term proportional to H . 
Without loss of generality, in this section we only consider interactions that include 
terms proportional to o x . These assumptions are often justifiable on experimental 
grounds |2*2^ l2*ïï| l2*4~] . The Hamiltonian of the system then becomes, 

H(t) = H + V(t) = -—°z + V(t)a x , (3) 



and the probability amplitudes evolve according to 

.d_ ( ai (t)\ = í-AE/2 V(t) \ ( ai (t) 
X dí \a 2 (t) ) { V(t) AE/2 ) \ a 2 (t) 



(4) 

Formal solution to (0J) may be written in terms of the time evolution operator U (t) 



as 



(:!?)-^(3)- 

In general, solving (J3J) and (jüj requires use of numerical methods. 
The time evolution operator Ü(t) may be expressed here as 

Ü(t) = Te-'Io 6 ^' = Te~ i Iol·^+ v ^) dt ' (6) 

00 (—\) n rt A rt A rt „ 

= T E^ / H(t n )dt n ... / H(t 2 )dt 2 / H^dh . 
n\ Jo Jo Jo 

The only non-trivial time dependence in U(t) arises from time dependent H(t) and time 
ordering T. The Dyson time ordering operator T specifies that H(ti)H(tj) is properly 
ordered: 

TH(ti)H(tj) = H(ti)H(tj) + 9(tj - U) [H(tj), H{t t ) . 

Time ordering imposes a connection between the effects of H(ti) and H(tj) and leads to 
observable time ordering effects [221123 121] Since time ordering effects can be defined 
as the difference between a result with time ordering and the corresponding result in the 
limit of no time ordering, it is useful to specify carefully the limit without time ordering. 
Removing time ordering corresponds to replacing T — ► 1 in @. 
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2.2. Analytical solutions 

In this paper we emphasize the utility of having analytic solutions, i.e. as comparec! 
to less transparent solutions obtained numerically. Unfortunately there are a limited 
number of conditions under which analytic solutions may be obtained. In order of 
increasing complexity these include: 

(i) Perturbative interactions [2H|. Here the interaction V(t) is sufficiently weak that 
the system largely remains in its initial state. The solution of (Jlj is trivial: 

U W~ {-ije- iAE ^V(t')dt' e" 

The mathematical validity condition is that the action associated with the external 
field is small, i.e. / V(t)dt « 1. 

(ii) Degenerate basis states |2*ïï] . In this case the energy levels of the two unperturbed 
states are nearly the same. For two-state systems the occupation probabilities are 
typically cos 2 (J V(t)dt) and sin 2 (J V(t)dt). Remarkably this form holds for both 
slowly and rapidly changing fields. Validity requires that the action difference 
associated with free propagation of the two unperturbed states be small, i.e. 
AEt « 1. 

(iii) Constant external fields. Here the interaction V(t) is a constant. The analytic 
solution, found from that for slowly changing fields given immediately below, is 
mathematically similar to the physically distinct RWA solution fZT\ . 

(iv) Slowly changing (adiabatic) fields j^H^H]- The analytic solution of (Jlj is, 



cos9(í) +iA|sine(í) -2i^sine(í) 
. -2ig|sin0(í) cos9(t)-iA| s ine(t) 

rt . 



Ú{t): 

where ü(t) = sj (AE) 2 + AV(t) 2 and 6(t) = £ Ü(t')dt'/2. The validity condition 
V(t) << íl 2 (t) is sometimes difficult to achieve. 

(v) RWA solutions. V(t) oscillates with a frequency u close to the resonant frequency 
of the transition between the basis states, ujq = AE. The RWA expression for U 
is the same as the expression for slowly changing fields given immediately above, 
except that Q(t) = üt, where ü 2 = V 2 + (Au;) 2 . The RWA is vàlid [27 when 
the frequency of the external field, lü, is nearly the same as that of the transition 
frequency, c^o, i.e. Au = u — ojq « ojq. 

(vi) A sudden pulse [2H] or series of sudden pulses (single or múltiple kicks). The 
bàsic validity condition [301 is that the external field is sharply tuned in time, 
i.e. AEt « 1, where r is the width of the pulse. This condition is met when 
t is relatively small. We examine this in detail below. In many experiment ally 
accessible cases one can build an external field using a combination of kicks. 

These solutions represent different scenarios, some of which can lead to a significant 
or even complete transfer of the population between basis states of the qubit. The 



Sudden switching in two-state systems 



6 



scenarios that allow complete transfer of the population are especially interesting for 
possible applications in the field of quantum information. Any measurement of a 
superposition state of the qubit leads to the collapse of the wavefunction and results 
in finding a qubit in one of its basis states. Therefore, the only states of a qubit for 
which one can predict the outcome of the measurement are the basis states themselves. 
Consequently, one needs a reliable way to drive the system into one of these states. Also, 
to form an arbitrary superposition state, one should be able to transfer any fraction of 
the population of the system into any of the states. This cannot be achieved with the 
perturbative scenario, for which the transfer of population is always incomplete, and 
some superposition states can never be formed. 

In the stationary or adiabatically changing field scenario, completeness of the 
transfer is limited by the ratio AE/V. Transfer is incomplete unless the energy levels are 
degenerate. In RWA, the probability amplitudes oscillate [2Zj with the Rabi frequency Cl, 
and completeness of the transfer can be adjusted by changing the detuning parameter, 
Au. Again, only in the limit of exact resonance, Auj — ► 0, is the transfer complete. 
Another technique, which is based on RWA and, when applicable, enables one to achieve 
almost complete transfer of the population, is STIRAP (stimulated Raman adiabatic 
passage [HI]). In this approach, a counter intuitive sequence of a pump pulse and a 
Stokes pulse is used to transfer the population via an intermediate state without losing 
any population due to the spontaneous decay of that state. This technique has proven 
to be very effective in a number of systems. The limiting factors there, however, are: 
i) it cannot be applied to degenerate or nearly degenerate systems, and ii) the pulses 
have to be applied adiabatically, which prevents one from using fast and ultrafast pulses 
(typical duration of pulses used in STIRAP is of the order of a nanosecond, whereas in 
the kick approach, the only restriction comes from the structure of the energy levels of 
the system, so picosecond and, in some cases, even femtosecond pulses can be used). 
In degenerate qubits, even higher degree of controllability can be achieved [2E!- And in 
some systems, a natural way to achieve a complete transfer of the population in a qubit 
is to apply a sudden pulse, or a kick. The focus of this paper therefore is on suddenly 
changing pulses, i.e. kicks, where population transfer can in some cases be complete. 



2.3. Pulses 

Kicks are an ideal limit of finite pulses, or sequence of pulses, each of some finite duration 
r. Each kick causes sudden changes in the populations of the two states. It is instructive 
to define phase angles for each individual pulse, namely, 

a = J V{t)dt , 

(3 = tAE/2 . (7) 

The angle a is a measure of the strength of the interaction V(t) over the duration of 
a given pulse. The angle j3 is a measure of the influence of H during the interaction 
interval r. 
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Exact analytical solutions can be obtained in the limit of kicks when the pulse 
applied at the time t — % becomes a 5-function, V(t) — > a^d^—t^), since the integration 
over time in (jHJ) becomes trivial. For a more realistic case where the pulse has a finite 
width, one may, at best, only obtain an approximate solution. 



3. Sudden switching 

In this section we consider single and múltiple kicks, where V(t) may be described in 
terms of delta functions in time, i.e. V(t) = J2k=i a k5{t — £&). We work primarily in 
the interaction representation, since the solutions are relatively simple and there are 
generally advantages with convergence in the interaction representation f22J. In the 
interaction representation the evolution operator has the general form 

Ü I (t)=Te- i tiWW 

= Texp(-ij\- i ^ Lt '^V(t')a x é^ t '^dt^ . (8) 
It's straightforward (e.g. using power series expansions) to show that 

g-i^í^gi^toz _ cos (/\Et)a x + sm(AEt)a y . 
Introducing a unit vector n(t) = {cos(AEt); sm(AEt); 0}, (jHJ) can be written as 

Ui(t) = Texp(-i^ V(t')n(t')-adA , (9) 

where a = {a x ; a y ; a z }. 

As mentioned above there are relatively few cases in which analytic solutions are 
available. Considered next is one set of such cases, namely singly and multiply kicked 
qubits. 



3.1. A single kick 

The bàsic building block is a two-state system subject to a single kick at time 
corresponding to V(t) = ctkS(t — í&). The integration over time is trivial and the time 
evolution operator in (Jüj) becomes 

s r . . / cosat — i sin a fc e~ lAEífe \ 

Uf (t = exp [-ia k n {t k )-a = .. h iAEtk 10 

y —i sm aike h cosafe J 

for t > tk- Here we used the identity e 1 ^ """ = cos0 / + isin0 o ■ u, where ü is an 
arbitrary unit vector. Note that Uf(t) is independent of t since the e lAEt factors, due 
to free propagation, are transferred from the evolution operator to the wavefunction in 
the interaction representation. 

Another way to evaluate Üj(t) is to use Ug(t) from the Schròdinger representation 
and to use the general relation, Uj(t) = e lHot Ug(t). For a single kick Us(t) has been 
previously evaluated [30J, namely, 

U SW ~ [ _ ie -iAE(t/2-t k ) sinak e -iAEt/2 CQSak ■ 
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There is an explícit dependence on time in Üg(t). Even in this elementary example the 
expression for the time evolution matrix, U, is simpler in the interaction representation 
than in the Schròdinger representation. 



For a kicked qubit initially found in the state ^ q J > the occupation probabilities 
are simply, 




P 1 (í) = |a 1 (t)| 2 = |t/ 1 fc 1 (í)| 2 = cos 2 a fc , 

P 2 (t) = \a 2 (t)\ 2 = \U^(t)\ 2 = sm 2 a k . (11) 



When the pulse width is finite, the corrections to (jllj) are 0{f3). These corrections 
result from the commutator of the free Hamiltonian Hq with the interaction V during 
the time r when the pulse is active. For example, in the case of a rectangular pulse of 
width r, the exact time evolution is given [30] by 



where a' = \/a 2 + fi' 2 . To leading order in (3, i.e. in the width of the pulse, the error in 
the kicked approximation is given by 



In the Schròdinger picture, time ordering effects are present even for a single ideal kick, 
specifically the time ordering between the interaction and the free evolution preceding 
and following the kick |3Uj . The time ordering effect vanishes in either the degenerate 
limit AEt — > or in the perturbative limit a — ► 0. 

In the interaction picture, time ordering effects disappear for a single ideal kick. 
This is easily understood by considering that in the interaction picture, time ordering 
is only between interactions at different times, Vj(í') and Vj(í"), not between the 
interaction V(t') and the free Hamiltonian Ho(t"), as in the Schròdinger case. For a 
single ideal kick, all the interaction occurs at one instant, and no ordering is needed. Of 
course, for a finite-width pulse, i.e. j3 ^ 0, time-ordering effects do begin to appear even 
in the interaction picture [301 ■ We note that the time ordering effect in the interaction 
picture is independent of the measurement time t, though it does depend on the pulse 
width t through the (3 parameter. 

3.2. Múltiple kicks 

Consider now a series of kicks, «i,a 2 , applied at t = ti,t 2 , ■■■,t n with t\ < t 2 < 

... < t n , i.e. a potential of the form V(t) = a\5(t — t\) + cv2S(t — £ 2 ) + ••• + o: n 5(t — t n ). 
In the interaction representation, this potential has the form: 




(12) 




(13) 



Vr(t) = {a x 5{t - t x ) + a 2 5(t - í 2 ) + ... + aj{t - t n ))n{t)-a . 
The evolution operator (JÜJ) becomes 




(14) 
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For a given order of pulses, the (|14|) can be written as 

jjmk _ eX p[—[ an fi(t n y ff] x ... x exp[— iain(ti)-a] , (15) 

which is a simple product of time evolution operators for single kicks. Using (fTUJ) . one 
obtains: 

ÜT k ={ COSan AEt - Íf * na » eriAEtn ) x ... (16) 
^-ísina^ 1 " cosa n J 

( cosai — i sin aie~ lAEtl \ 

\ —i sin aie lAEtl cos ot\ ) 

This can be evaluated for an arbitrary combination of kicks, so the analytical expression 
for the final occupational probabilities for the basis states can be obtained. As we shall 
explicitly demonstrate later, the order in which the kicks occur can make an observable 
difference. 

3.2.1. Two arbitrary kicks The simplest example of a series of arbitrary kicks is a 
sequence of two kicks, of strengths «i and a 2 , applied at times t x and í 2 - Then (fTüj) is 
easily solved, namely, 

Üf ] = Uf x Üj 1 =e^[-ia 2 n{t 2 )a}e^[-iain{t 1 )a} 



f X 


u^ = 






u 21 


TJ* 



;i7) 



where 



-i/\Et. 



Un = cos «i cos «2 — sin ai sin a 2 e 
U 2 i = — ie 2 + (cosaiSina 2 e 2 " +sm« 1 cosa 2 e 2 J • 

Here í_ = t 2 — ti, and t + — t\ + £2- In the limit t 2 — > íi, p7)l reduces to (frü|) with 
a — > «i + a 2 . Note that [Üf 2 , Üf 1 ] 7^ so that the time ordering of the interactions is 
important. In the interaction representation the expression for the matrix elements of U 
and the corresponding probability amplitudes are a little simpler than the corresponding, 
physically equivalent, expressions in the Schròdinger representation, which include an 
unnecessary explicit dependence on time. This reflects the idea that the interaction 
representation takes advantage of the known eigensolutions of H Q . 

The àlgebra for doing a combination of two arbitrary kicks, one proportional 
to o y and the other proportional to a x , is very similar to that for two arbitrary 
kicks proportional to a x . For a single a y kick, similarly to (jlOj) one quickly finds 
jj ky _ ( cos a k - e ~ íAEt k si n ak \ 

\e lAEtk smak cosa^ I 

Then using JJ k2xkly = Ü k2x x Ü kly , one finds that the matrix elements for a a y kick 
at íi followed by a a x kick at t 2 are, 

Un = cosai cos a 2 — i sine*! sin a 2 e~ lAEt ~ , (19) 
U 2 i = e 2 + (cosa2 sinaie 2 ~ — i sin a 2 cos aie 2 ~) . 



Sudden switching in two-state systems 



10 



From that, the transition probability to the off state is, 

p k 2x k ly = p 2i \2 (20) 

= cos 2 «i sin 2 a 2 + sin 2 oc\ cos 2 a 2 + - sin 2a x sin 2a 2 sin AEt_ . 

If the a x and er^ kicks are reversed, then í_ — > — í_ so that pf xkly — p£ 2yklx = 
sin 2a;i sin 2a 2 sin AEt_ . This difference oscillates between ±1 when a x = «2 = 7r/4. 
This effect can be observed. In contrast there is no observable difference for two o~ x 
kicks as may be easily shown from (|18|). 



3.2.2. Two identical kicks If the pulses for two kicks are identical (ai = a2 = a), 
then ()18j) simplifies further to 

U u = cos 2 a - e- iAi?í - sin 2 a , (21) 

,AE, AE 

Un = — le 2 r + sin 2a cos -— í_ . 

In the limit í 2 — * h ÍÍ2*T|) reduces to (fTÜj) with a doubled. 

Similarly, for two a x kicks of equal magnitude but opposite sign (ai = — a2 = a) 
applied at times t\ and t% y 

U u = cos 2 a + e~ iAEt - sin 2 a , (22) 



■AE, AE 

U21 — — smzae 2 + sm — — í_ 



2 

which in the limit Í2 - ► íi reduces to the identity matrix. 



3.2.3. Three arbitrary kicks For three arbitrary a x kicks of strengths ai, a2, and a^, 
applied at times ti, Í2, and Í3, 



Üf = Uf x Uf x Uf 
(Un -U* 21 \ 
U« U* u ) ' 



where 



(23) 



Un = cos ai cos a2 cos — sin ai sin a2 sin (24) 
x ( e iAE ^-^ cot ai + e iAE ^~^ cot a 2 + e iAi ^- Í2 > cot a 3 ) , 

U 2 i = i(sin ai sin a 2 sin a 3 e lAE ( Í3_Í2+íl ) — cos ai cos a 2 cos a 3 

x (e iAEíl tan a x + é AEt2 tan a 2 + é AEt3 tan a 3 )) . 



4. Calculations 

In this section we present the results of numerical calculations of occupation probabilities 
for transitions in a two-state system, caused by a series of narrow Gaussian pulses of 
width r. We study the effects of the order in which pulses are applied on the final 
occupation probabilities. First we illustrate the results obtained in the previous section 
using a model two-state system. Then we present realistic calculations for 2s — > 2p 
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transitions in atòmic hydrogen. We also discuss in detail the applicability of a two-state 
model to this transition. 

4-1. A model two-state system 

Here we present the results of numerical calculations for transitions in a model two-state 
system. We directly integrate flU) using a Standard fourth order Runge-Kutta method. 
In our calculations we use an interaction of the form Vk(t) = (a/ c /y7r'r)e~ ( - — ^ , i.e. we 
replace an ideal kick (delta function) by a Gaussian pulse of strength a k centered at t k 
with width r. When r is small enough for the sudden, kicked approximation to hold, 
this should give the same results as the analytic expression for a kick above, i.e. as 
r -.0, V k (t)^a k 5(t-t k ). 

4-1.1. Two similar kicks First we calculate the probabilities for transitions caused by 
two similar kicks. Both pulses are proportional to o~ x , but the action integral vàlues are 
different: V(t) = Vi(t)-a x +V2(t)·a x (cf. (|I7J)). In figure[T]we show results of a calculation 



very nearly achieved since we choose r to be a factor of 10 3 times smaller than the Rabi 
time, Ta, for the population to oscillate between the two states. Peaks and dips in the 
Pzit) graphs, that occur during application of a second pulse, reflect the following fact. 
By the time second pulse is applied, the system already is in a superposition of the 
bàsic states, and it's this superposition that undergoes a precession when the pulse is 
on. The final occupation probability of the target state doesn't depend on the order in 
which the kicks are applied. In figure El the results of a similar calculation for broader 
pulses of width r = 0.0057a are shown. The shape of the P(t) graph reflects the fact 
that the pulses have finite width. However, the outcome of the calculation still doesn't 
depend on the order in which the kicks were applied. 

4-1.2. Ordering effects As we have shown in the previous subsection, if two pulses act 
on the system, then the outcome of the process is independent on the order in which they 
are applied as long as the pulses are similar (e.g. two a x or two o~ y pulses). However, if 
the two pulses that act on the system have different structure (e.g., one is a o x pulse, 
and the other one is a a y pulse), then the results can be significantly different if different 
sequences of pulses are used. In figure the occupation probability of the target state is 
calculated using two different sequences: an a.\0~ x pulse followed by an a.io y pulse (sòlid 
line) and an a^Oy pulse followed by an a.\0~ x pulse (dashed line). All the parameters are 
identical to the ones used in the previous part (two a x pulses), except for the structure 
of the pulses. The difference between two sequences is obvious. The effect of using 
different order of pulses, as well as the occupation probabilities for each case, are in 
very good agreement with the vàlues calculated analytically (cf. püjl ). 



for the probability P^if) that a system initially in the state 




makes a transition into 



the state 




when perturbed by two pulses applied at t% and £2. The ideal kick is 
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Figure 1. Target state probability as a function of time. Two <j x kicks applied at 
t\ and t 2 . Width of the kicks: r = O.OOITa- Action integral vàlues: a.\ — 0.17T, 
a 2 = 0.157T (chosen arbitrarily). The sòlid line: probability for a\ followed by a 2 ; the 
dashed line: probability for a 2 followed by a±. The sharp dip due to using a pulse of 
finite width is explained in the text. The final probability doesn't depend on the order 
of kicks. 
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Figure 2. Target state probability as a function of time. Same as figure ^ but 
broader kicks (kick width r = 0.005Ta). The sòlid line: probability for ct\ followed 
by a 2 ; the dashed line: probability for a 2 followed by a.\. The final probability still 
doesn't depend on the order of kicks. 
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Figure 3. Target state probability as a function of time. Two kicks, u\o x and a2<? y , 
applied at t\ and t 2 - Kick width r = O.OOITa- Action integral vàlues: otx = 0.17T, 
012 = 0.157T. The sòlid line: probability for ot\o x followed by oi2(Jy\ the dashed line: 
probability for 0.1a y followed by ot\<j x . The final probability depends on the order in 
which the kicks are applied. 



Even for a sequence of pulses of the same structure, the order of pulses can be 
significant. To illustrate this fact, consider a series of three a x pulses. The results of 
numerical calculations are shown in figure |U As long as the time intervals between 
pulses are not the same, the outcome of the process does depend on the sequence in 
which the pulses are applied. 

4-2. 2s — 2p transition in hydrogen 

In this section we present the results of numerical calculations for 2s — > 2p transitions in 
atòmic hydrogen caused by a series of Gaussian pulses of width r. Applicability of a two- 
state approximation to this transition is discussed in detail in the Appendix. Specifically, 
we consider the fine structure splitting of the 2p state (target state) into 2pi/ 2 and 2p 3 / 2 
states. As we show, the phase difference accumulated during free evolution between 
pulses oscillates with the period T r = 2ir/Ef s , where Ef s ~ 10956 MHz is the fine 
structure splitting. Therefore, the same superposition state is formed periodically. By 
choosing time intervals between kicks to be integer múltiples of T r , one can effectively 
treat the superposition of 2pi/ 2 and 2p 3 / 2 states as one state (2p state), that is coupled 
to the initial 2s state. The occupation probabilities of the initial state 2s and the target 
state 2p are evaluated by integrating two-state equations using a Standard fourth order 
Runge-Kutta method. This enables us to verify the validity of our analytic solutions for 
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Figure 4. Target state probability as a function of time. Three a x kicks applied at ti, 
t 2 , and t 3 . Kick width t = O.OOITa- Action integral vàlues: a% = O.Itt, a 2 = 0.157T, 
Q!3 = 0.257T. The sòlid line: probability for u\ followed by a 2 followed by the 
dashed line: probability for «3 followed by a 2 followed by ai. The final probability 
depends on the order of kicks. 



kicked qubits in the limit as r — > and also to consider the effects of time ordering. 

The splitting between the 2s and 2pi/2 states in atòmic hydrogen (Lamb shift) is 
AE « 1057 MHz. The corresponding time scale (the Rabi time that gives the period of 
oscillation between the states) is Ta ~ 10~ 9 seconds. This gives rise to the first limitation 
on the duration of the pulse, r: it has to be significantly smaller than Ta, otherwise the 
pulse will not be sudden and the kicked approximation will fail. On the other hand if r 
is too small, then the interaction will have frequency components that couple the initial 
state to other states. Specifically, if 1/r is greater than (E^ p — E2 S ) ~ 10 15 Hz, then the 
interaction will induce transitions into states with n > 3 and the system will not be well 
approximated by a two-state model. Also there is another constraint in our case. If the 
time of interaction becomes comparable to the lifetime of one of the active states (the 
less stable 2p state has a lifetime of ~ 1.6ns), then the dissipation effects (spontaneous 
decay into the lower states outside two-state model) cannot be neglected. 

Here we use Gaussian pulses with width r = lps, and limit the time of interaction 
to ~ 600ps. Single and múltiple pulses of such width (and even much shorter) are 
achievable experiment ally (e.g. half-cycle electromagnètic pulses, [33 ). We explicitly 
include in the numerical calculations both 2pi/2 and 2^3/2 states, and keep the time 
separation between the pulses t 2 — t\ = T r . The loss of population from the two-state 
system due to spontaneous decay (2p — > ls) is also included. The evaluation of a in 
terms of the dipole matrix element for the 2s — 2p transition is discussed in a previous 
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Figure 5. Target state probability as a function of time for 2s — > 2p transition 
in atòmic hydrogen. Two a x kicks applied at t\ — 20ps and t 2 = 593. 5ps. Width 
of the kicks: r = lps. Action integral vàlues: u\ — O.lvr, a 2 = 0.15tt. The sòlid 
line: probability for ol\ followed by a 2 ; the dashed line: probability for a 2 followed by 
a\. The only difference in the final probability for different orders of kicks is due to 
dissipation. As in figure the sharp dips are real and due to the finite width of the 
pulsc. 



paper j2Ej- We present results for the occupation probability of the target state, P 2 , 
which includes both 2pi/ 2 and 2p 3 / 2 states, as a function of time. 

4-2.1. Two similar kicks In figure 03 we show the results of a calculation for the 
probability P%(t) that a hydrogen atom initially in the 2s state makes a transition 
into the 2p state when perturbed by two a x Gaussian pulses applied at t\ and t 2 . We 
have obtained our results by numerically integrating the coupled equations, 
/oi\ / AE -V(t) -V2V(t)\ (a x 

\a 3 J \-y/2V(t) {E fs -ïl)) \a 3 , 
where T w 626 MHz is the decay rate for the 2p state, and we set E 2pi/2 = 0, so 
-^2p 1/2 = AE and E 2v , il2 = Ef s . For the calculation, we used the following parameters: 
kicks applied at t± = 20ps and t 2 = 593. òps (separation between the pulses is equal 
to the revival time T r for the superposition of 2pi/ 2 and 2p 3 / 2 states); action integral 
vàlues: a,\ = O.lir, a 2 = 0.157T. The final occupation probability of the target state 
doesn't depend on the order in which the kicks are applied. The small difference in the 
final probabilities is entirely due to dissipation effects since the decay rates of 2s and 
2p states in hydrogen differ by nine orders of magnitude. Removing dissipation yields 



«2 = ~V(t) -if • a 2 , (25) 
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Figure 6. Target state probability as a function of time. Two kicks, a\a x and a2<? y , 
applied at ti = 20ps and t 2 — 593. bps. Kick width r = lps. Action integral vàlues: 
a\ = 0.17T, a.2 = 0.15tt. The sòlid line: probability for a\a x followed by a.20 y \ the 
dashed line: probability for ot2<J y followed by ot\a x . The final probability depends on 
the order in which the kicks are applied. 



results indistinguishable from figure 

In figure |ül the occupation probability of the target state is calculated using two 
different sequences: an a.\<7 x pulse followed by an a2& y pulse (sòlid line) and an a.iG y 
pulse followed by an a.\O x pulse (dashed line). All the parameters are identical to the 
ones used in the previous part (two a x pulses), except for the structure of the pulses. 
Now the difference between the effects of two sequences is obvious. 

4-3. Effect of time ordering in a doubly kicked system 

In this subsection we use our analytic expressions to examine the effect of the Dyson 
time ordering operator, T, in a kicked two-state system. The effect of time ordering 
has been considered previously in the context of atòmic col·lisions with charged partides 
|18| l%m El E2 I2H1 E3] arid differs somewhat from the order in which external pulses 
are applied, as illustrated below. As is intuitively evident, there is no time ordering in 
a singly kicked qubit |30j since there is only one kick. The simplest kicked two-state 
system that shows an effect due to time ordering is the qubit kicked by two equal and 
opposite pulses separated by a time í_ = t 2 — t\. The evolution matrix Uj kk for this 
system of may be rewritten for convenience (as may be easily verified) as, 
e'' 1 ^ 1 - (cos ^-t- + i cos 2a sin 4f t_) e -^Et+ gin 2a s j n 

_ e iAEt+ gin 2a S i n e 1 ^*- (cos - i cos 2a sin 4f £_) 
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The limit of no time ordering, i.e. T — > 1 in may in principle be generally 
obtained [3D] by replacing J^V^dt with Vt, where V is an average (constant) value 
of the interaction. In the case of two kicks of the same magnitude and opposite signs, 
it is then straightforward to show that, 

fM-kk _ -m _ ( cos(2« sin e ~ iAEt + sin(2« sin \ 

7 l> ~~ V-e iAEt +sin(2asin^í_) cos(2« sin fí_) 1 J 

In this example we now have analytic expressions for the matrix elements of both Uj 
that contains time ordering and Üj°^ kk that does not include time ordering. 

Let us now pause to examine the difference between time ordering and time reversal 
in this simple, illustrative example. Reversal of time ordering means that, since a^ and tk 
are both reversed, both t_ — > — í_ and a — > —ct. In this case one sees from the equations 
above that Üj kk is not invariant, i.e. phase changes occur, while Üf^ kk remains the 
same. Thus Üf kk changes when the time ordering is changed, but Üf^ kk does not 
change. For time reversal p5] t± — > —t± and, since the initial and final states are also 
interchanged, U W. In this case one sees by inspection of the above equations that 
both Üf kk and Üf ] ~ k ' are invariant under time reversal, as expected. As shown below 
when the symmetry of the kicks is broken the difference between Ü k Ü k ' and Ü k 'Ü k can 
be observed. 

We have shown above both algebraically and numerically that for two kicks 
proportional to a x , the order of the kicks does not change the final population transfer 
probability Pi. However, interestingly, this does not mean that there is no effect due 
to time ordering in this case. As we show next, there is an effect due to time ordering 
in this case, even though reversing the order of the kicks has no effect. The effect 
of time ordering on the occupation probabilities may be examined by considering the 
probability of transfer of population from the on state to the off state with and without 
time ordering, namely, 



P 2 = \U 12 \ 2 = |sin2a sin^-í_| 2 = |esin0| 2 , (28) 

pp = \U$\ 2 = |sin(2asin^í_)| 2 = |sine0| 2 , 

where e = sin 4pt- and <fi = 2a. 

The effect of time ordering is shown in figure where P 2 — P^ ^ is plotted as a 
function of = 2a, corresponding to the strength of the kicks, and e = sin 4pí_, which 
varies with the time separation of the two kicks. The effect of time ordering disappears in 
the example we present here in the limit that either the interaction strength or the time 
separation of the pulses goes to zero. For small, but finite, vàlues of both the interaction 
strength and the time separation of the pulses, the effect of time ordering is to reduce 
the probability of transition from the initially occupied state to an initially unoccupied 
state. That is, in this regime time ordering reduces the maximum transfer of population 
from one state to another. As either of these two parameters gets sufficiently large, the 
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Figure 7. Diffcrence in population transfer probability, P2 — P 2 vs - e = sm ^r^- 
and (j> — 2a. Here í_ = t 2 — t\ is the time between the pulses and a = j V(t')dt' is a 
measure of the interaction strength. The two-state system is kicked by a sharp pulse 
of strength a at time t\ and by an equal and opposite pulse at time Í2. The difference, 
P2 — í*2 , ^ s ^ue t° time ordering in this qubit. 



effect of time ordering oscillates with increasing vàlues of the interaction strength or the 
separation time between the two pulses. Time ordering effects are present even though 

jj—kk jjk—k 

5. Discussion 

Clearly one may extend this approach past two kicks or three. Since arbitrary pulses 
can be built from a series of kicks, in principle one may build arbitrarily complex pulses 
using kicks. While adding more kicks is straightforward, the àlgebra becomes more 
difficult. Also, the number of natural systems in which the validity conditions apply 
diminishes as the pulse becomes more complex. Hence it may be sensible to seek cases 
that have sufficient symmetry so the analysis is both simple and applicable. It has been 
previously noted (HO], for example, that a simple expression, corresponding to Floquet 
states, exists for a periòdic series of narrow pulses. 

Part of the motivation for this paper grew out of an effort to define correlation 
in time [211], based on effects of time ordering. While in principle we have found no 
fundamental problem with this effort, we have found that it is often difficult to find 
expressions for the time evolution matrix in the limit of no time ordering, namely U°. 
Furthermore, U° can depend on the representation used [3Ü]. We also note that time 
ordering may occur in the degenerate limit, e.g. when H = H cr z + Vi(t)a x + V^i)^. In 
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recent applications in collision dynamics using perturbation theory |2T| . time ordering 
is removed by use of degenerate states since the external interactions do not contain 
more than one type of coupling. But in these calculations there is no difference between 
time ordering and time correlation. The most reliable way to remove time ordering 
is replacement [30] of the instantaneous interaction V(t) by its time averaged value, 



In summary, analytic solutions for two-state systems (e.g. qubits) strongly 
perturbed by a series of rapidly changing pulses, called 'kicks', have been developed 
and discussed. Such analytic solutions provi de useful physical insight, which together 
with more complete numerical methods [3U ES] may be used to solve more complex 
problems. For a series of kicks the evolution matrix may be expressed as a time ordered 
product of single kicks. We have explicitly considered in detail single, double, and triple 
kicks. While there is no difference in the population transition probability if two o x kicks 
are interchanged, time ordering does have an observable effect. The effect happens to 
be the same for both of these orderings. If a a x kick is interchanged with a a y kick in 
a doubly kicked system, the difference can be observed in most cases. If three a x kicks 
are used, different orderings can also be observably different. The effect of using pulses 
of finite widths has been studied numerically for 2s — 2p transitions in atòmic hydrogen. 
If the pulse width is much smaller than the Rabi time of the active states, then the 
analytic kicked solutions are vàlid. Such pulses can be created experimentally using 
existing microwave sources. The difference between time ordering and time reversal 
has been specified. Time ordering can have observable effects. Under time reversal 
the quantum amplitudes are generally invariant. Our results may be extended to an 
arbitrary number of kicks. However, without simplifying symmetries, the solutions 
become more complex, and the applicability of this approach becomes more limited, as 
the number of kicks increases. 

Appendix 

Here we discuss some details concerning the use of the two-state approach for studying 
the 2s - 2p transition in hydrogen. First, we note that the eight n = 2 states (two 2sy 2 
states, two 2p\/ 2 states, and four 2p 3 / 2 states) are nearly degenerate and well separated 
in energy from states with n ^ 2. Thus, smooth external pulses may easily be chosen 
long enough to prevent field-induced population transfer out of the n = 2 subspace, and 
requiring us only to include the spontaneous decay rate T from 2p to ls. Furthermore, 
rotational invariance around the axis of the external elèctric field E leads to conservation 
of total angular momentum component along that direction, allowing a given 2s\/ 2 state 
to couple only to one 2py 2 state and one 2p 3 / 2 state. Specifically, starting with an initial 
2s\/ 2 state with spin polarization at some angle % relative to E, the accessible subspace 
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is spanned by the three basis vectors 

|2 S ) =\£ = 0, m = Q)\ X ) 
\2p) =\£=l,m = 0)\ X ) 
|2j/) = cos||i = l 



(A.l) 



m 



+l)|I)+sin||£ 



1 , m 



-1)1 T) 



where the initial spin state is \x) = cos|| |) +sin^| j.}, and both orbital angular 
momentum and spin components are measured along the direction of E. 

Since the external pulse does not change the orbital angular momentum component 
m, the external field couples only the two states \2s) and \2p) in the above basis, e.g. 

/ 1 \ 

(A.2) 



V(t) = V(t) 












for a a x pulse. The free Hamiltonian H is diagonal in the basis of good total angular 
momentum j. In the basis of (|A.1J) . 
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2 rp 

3- C/2 P3/2 



V o 





3 a f s 



\E 







\ 



2P3/2 



+ ÍE, 



(A.3) 



2P1/2 



where Ef s — E 2p3/2 — E 2pi/2 is the fine structure splitting, and we take the decay rate T 

to be the same for 2py 2 an d 2p 3 / 2 . 

For narrow pulses, whose inverse width 1/r is large compared both with the splitting 

AE between the 2si/ 2 and 2p\/ 2 energies (Lamb shift) and also Ef s (fine structure), the 

free propagation may be neglected during the time of the pulse. Then the full propagator 

in the interaction representat ion may be written as a product of kick operators of the 
form e +ítH t ne -iJdtv n (t) e -iH t„^ where tn ig the time of the n th kick Now e -i/díV„(í) ig 

block-diagonal by construction ()A.2|) . with the 2p' state decoupled. In between pulses, 
amplitude oscillates between the 2p and 2p' states with period T r = 2n/Ef s . However, 
if we now choose all inter-pulse spacings to be integer múltiples of this period, 



Aí n t n j r \ t n TfiT r , 



then the free propagation between kicks also becomes diagonal: 



-iH At n 



\ 



e -iE 2si/2 At„ 








(-i£ 2 p 1/2 -r/2)Aí n 








e (- Í ^Pl/2- 



■r/2)At„ 



(A.4) 



(A.5) 



as may easily be checked explicitly by writing the above free propagator in the 2pi/ 2 , 
2p 3 / 2 basis and noting that the two basis vectors acquire the same phase e lE2 Pi/2 Atn = 
e lE2p 3/2 tn . Thus the 2p' state decouples entirely and its occupation probability will 
always be zero when we view the dynamics stroboscopically with period T r starting with 
the time t\ of the first kick. The three-state dynamics therefore reduces to two-state 
dynamics in the 2s, 2p subspace. 
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Finally, as long the the inter-kick spacings are all integer múltiples of T r , we may 
also use the two-state formulas presented in the main body of the paper to evaluate the 
occupation probabilities at an arbitrary time between kicks or after the last kick. We 
need only remember that the 2p probability that we compute at these arbitrary times 
is the total probability for being in either the 2p and 2p' state, or equivalently the total 
probability for being in either the 2py 2 or 2p 3 / 2 state. 
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